Experimental synchronization of spatiotemporal chaos 
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We report the first experimental evidence of a successful synchronization of spatiotemporal chaos. 
The experiments were performed on two unidirectionally coupled, nonlinear-optical systems of single- 
feedback type. The synchronization was investigated for different degrees of complexity of the 
spontaneous structures. In all cases, the cross-correlation between the two system states was found 
to increase with the strength of the coupling. Numerical simulations yield comparable results and 
throw a light on the role of spatial inhomogeneities, which hamper a perfect synchronization. 

PACS numbers: 89.75.Kd, 05.45.Xt, 47.54.-|-r, 42.65. Sf 



It is well-known that two coupled chaotic oscillators 
can synchronize P, Q] . The extension of this concept to 
spatially extended systems is a rather new topic. Sponta- 
neously forming complex and dynamic spatial structures 
("spatiotemporal chaos") play a role in numerous fields, 
from biology over chemistry to different areas of physics 
0, e.g. fluid dynamics, solid state physics, or nonlinear 
optics While as bridge to the temporal phenomena, 
often spatially discrete system like coupled oscillators are 
regarded, we will here focus on continuous systems. 

To achieve synchronization between two separate sys- 
tems, the system states must be communicated between 
the both. In most pattern forming systems, such as in 
fluiddynamical or chemical systems, the experimental re- 
alization of this communication is rather difficult. The 
spatial distribution of physical quantities, like a flow field 
or the concentration of chemical compounds, must be 
manipulated dynamically by injecting an external sig- 
nal. Therefore, up to now only numerical simulations 
of coupled prototype models are performed, mostly in 
one spatial dimension only 0, S S Hi 01 • Both uni- or 
bidirectionally coupling, and the synchronization of not 
completely identical systems ^ are regarded. Following 
concepts of controlling-chaos, the coupling signal is some- 
times derived from a comparison of both system states, 
e.g. is the difference of two quantities. Since the transfer 
of the full spatial distribution of a system state is con- 
sidered to be a problem, often the possibility to use only 
few discrete coupling channels is investigated 0, 01 • In 
the present report instead, we regard the coupling of the 
full, space- and time-continuous state of spatially two- 
dimensional systems. 

Light waves are ideal as carrier for the coupling sig- 
nal, since they can carry almost arbitrary spatial and 
temporal profiles. Hence, it almost suggests itself to 
use nonlinear-optical pattern forming systems, where the 
light field is already a central physical quantity. Our 
experiment belongs to the class of the so-called single- 
feedback systems, which have become quite popular over 
the last years [ll IH El El El . 

In our system, a reflective Liquid Crystal Light Valve 
(LCLV) is used as optical nonlinearity, providing a self- 
defocusing, saturable Kerr-type nonlinearity. Its large 
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FIG. 1: Schematic experimental setup. For clarity, the pump 
beam is from the beginning already drawn as the two beams 
of the subsystems A and B. Components for beam expansion, 
for detection and for the bias beam of system A are left out 
for simplicity. More details are given in the text. 



nonlinear sensitivity allows us to realize large aspect ratio 
patterns with moderate laser powers. The LCLV consists 
of two thin layers, separated by a mirror and sandwiched 
between two transparent electrodes to which a low ac 
voltage is applied '15|. One layer is a photoconductor 
(L'C), which changes its conductivity according to the 
intensity profile l^{x,y) of an incident light wave. This 
results in a space dependence of the voltage drop over 
the second layer, a liquid crystal (LC), affecting its effec- 
tive refractive index. A light wave, passing the LC layer 
{read side) and being internally reflected, consequently 
acquires a transverse phase shift profile, which is deter- 
mined by the intensity profile at the PC {write side). 

The LCLV is put into an optical feedback loop: an 
almost plane pump wave is first phase-modulated and 
refiected by the LCLV read side. The modulated wave 
is then fed back to the intensity sensitive write side by 
means of beam splitters (BS), mirrors (M) and lenses 
(L), as indicated in Fig.^ During the free space propa- 
gation through the feedback loop, diffraction causes the 
phase profile to be transferred into an intensity modula- 
tion, such closing the feedback. For detection purposes, 
a fraction of the feedback wave is coupled out with the 
beam splitter BS2. We use a digital camera to record an 
intensity distribution Iw{x,y,t) equivalent to the one at 
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FIG. 2: Snapshots of spatiotemporally chaotic states without 
coupling (7 = 0) in inverse gray scale (dark corresponds to 
high intensity) at 1.5 times pattern forming threshold with 
closed low-pass filter (upper panel), and with open low-pass 
filter (lower panel). A and B denote master and slave system, 
respectively. 



the LCLV's write side. More details about this system 
and the theoretical modelhng can be found in jl^l . 

Above a certain threshold of the pump intensity, the 
uniform, plane-wave solution becomes modulationally 
unstable with respect to a critical transverse wave num- 
ber kc- At higher respective thresholds, also higher or- 
der critical wave numbers become unstable. Increasing 
the pump intensity from below threshold, first a sta- 
tionary hexagonal pattern develops 0, 0|, and then 
the pattern becomes increasingly disordered and dy- 
namic, finally turnin g in to a very complex spatiotem- 
porally chaotic state [IgI llTj . Typical snapshots of dy- 
namic experimental structures are shown in Fig. [2| Just 
above threshold, remainders of perfect order are visible 
as hexagonal domains, the typical size of which shrinks 
with pump intensity. 

Our system includes a spatial low-pass filter in the 
feedback (aperture Al in the Fourier plane between the 
lenses LI, L2). By setting the low-pass cutoff just above 
the first critical wave number, the evolution of spatial 
chaos is impeded and the amount of disorder grows more 
smoothly with the pump intensity 0,0. With an open 
low-pass filter instead, pronounced disorder sets in al- 
ready closely above the primary threshold. Hence, the 
closed low-pass allows us a finer adjustment of the de- 
gree of complexity. 



Because of their versatility, LCLV feedback systems are 
frequently used to investigate spontaneous optical struc- 
tures. In some configurations, the feedback wave is ro- 
tated or shifted laterally [l^, [3 . In the present version 
instead, the full rotational and translational symmetry is 
retained and diffraction represents the dominant spatial 
coupling mechanism. 

We have recently shown, how this system can be syn - 
chronized by injecting perfect hexagonal patterns |19l |. 
The question addressed in this report is, whether such 
a synchronization can also be realized for spatiotempo- 
ral chaos. For this purpose, we regard two systems, one 
of which runs autonomously (master A), while the other 
(slave B) is exposed to the attenuated signal generated 
by the master system. Since it is difficult to get two iden- 
tical LCLVs, we divided the active area of a single LCLV 
into two independent subsystems by inserting a mask A2 
with two circular holes (each diameter = 4 mm). Be- 
cause of inhomogeneities of the LCLV, the two systems 
A and B cannot be expected to be perfectly identical. 

In the feedback, a fraction of the light waves was ex- 
tracted with the beam splitter BS3. While the signal 
from the slave system B was blocked, the wave from the 
master A was injected into B with the beam splitter BS4. 
A halfwave-plate A/2 and a polarizer P were used to at- 
tenuate the injected signal, i.e. to set the coupling factor 
7. The perpendicular polarization behind the polarizer 
prevented interference effects. The recorded intensity dis- 
tributions do not contain this coupling signal, because it 
was injected behind beam splitter BS2. 

Due to the injected signal from A, the feedback inten- 
sity in B changes to I^b I'^b = ^wB + lIwA- With- 
out further modifications, system B would experience a 
total feedback power larger than system A. As conse- 
quence of the saturation of the nonlinearity, this would 
cause differences between systems A and B. In order 
to reduce this effect, a uniform bias beam was superim- 
posed onto A, such that I^a I'wA — ^wA + hias with 
hias — i{Iwa), where the brackets (•) denote spatial av- 
eraging. The seemingly simpler attenuation of I'^g was 
hindered by technical reasons. 

With the spatial low-pass filter cutoff set just above 
the critical wave number, the pump intensity Ip was set 
to different values above the pattern forming threshold 
Ith- An additional measurement series was carried out 
with open low-pass and the pump at Ip/Ith = 1-5. In 
each particular measurement, the coupling strength was 
set to a fixed value and a temporal image sequence of up 
to 100 s duration was recorded. From each image, we 
computed a cross-correlation function between A and B 

r,,A A ^\ {IwA{x,y,t) ■ i^six - Ax,y - Ay,t)) 
C{Ax,Ay,t) = , , 

(1) 

where I^a and I^b are the experimentally recorded 
intensity distributions of master and slave system, re- 
spectively. The tildes stand for the deviation from the 
mean value / = /—(/). In order to exclude bound- 
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FIG. 3: Time- averaged cross-correlation coefficient {cAB)t, 
plotted over the coupling strength 7. □: pump intensity set to 
1.5 times pattern forming threshold, O: to 3 times threshold, 
both for only one critical wave number active (closed low- 
pass); A: 1.5 times threshold with all critical wave numbers 
included (open low-pass). 



ary effects, only the central parts of the active areas 
(85% of the diameter) have been considered here. As 
measure for synchronization, the correlation coefficient 
CAsit) = C(0, 0,t) was used. Fig. |31 shows the time- 
averaged correlation coefficients, plotted over the cou- 
pUng factor 7. The errorbars indicate the RMS value 
of the temporal variation during the recorded sequence. 
We have checked that the cross-correlation coefficients do 
not show temporal drifts, which confirms that the cou- 
pled systems were in their asymptotic states. 

The increase of the cross-correlation coefficient cab 
with the coupling strength 7 clearly substantiates that 
the two systems become synchronized by the coupling. 
Because of the experimental imperfections and remain- 
ing differences between the systems A and B, the corre- 
lation coefficient never reaches unity. When the coupling 
parameter is larger than about 7 > 0.5, the correlation 
drops again. This behavior can be assigned to the fact 
that the systems become increasingly different, because 
the bias intensity for the master system A is an insuffi- 
cient way to balance the unidirectional coupling. Fig|31 
also shows that the cross-correlation drops with increas- 
ing distance from threshold, meaning that the synchro- 
nization is harder to achieve for stronger disorder. 

Results very similar to the cross-correlation were ob- 
tained by computingthe mutual information and the syn- 
chronization error pf, i.e. the normalized mean square 
difference between I^a and I^b- A snapshot of a syn- 
chronized state is presented in Fig.^ for the pump inten- 
sity just above threshold and closed low-pass filter. For 
weaker coupling, the correlations are smaller and difficult 
to be recognized by eye. 

Without low-pass filtering, a larger number of critical 
wave number bands contributes to the pattern formation 
process, resulting in very complex and dynamic struc- 



FIG. 4: Snapshot of synchronized states with coupling factor 
7 = 0.5 at 1.5 times pattern forming threshold (low-pass filter 
closed) . 



tures (cf. Fig. [3 lower panel). In this case, the max- 
imum achievable correlation is much lower and already 
saturates for 7 > 0.2 (cf. Fig. The reduced correla- 
tion is probably the consequence of the increased sensi- 
tivity to perturbations. Without low-pass filtering, the 
system supports a much broader bandwidth and conse- 
quently high-frequency spatial noise, such as speckles, 
gain more influence. However, the increase of correla- 
tion at lower coupling strengths still gives evidence for a 
partial, coupling-induced synchronization. 

In purely temporal systems, a so-called lag synchro- 
nization can be found, where the synchronized system 
follows the master system with a certain time lag 111, 13] . 
In a spatially extended system this would correspond to 
a lateral shift (or rotation?) between the two patterns. 
Such a shift would nota bene be connected to a symme- 
try breaking by the selection of a particular direction. 
We find, however, that the dominant peak of the cross- 
correlation function Eq. |^ does not change its location 
{Ax = Ay = 0). Hence, there is no lateral shift between 
the synchronized structures. 

For comparison, numerical simulations of the full 
model equations were carried out. Details on the nu- 
merical algorithm are given in pol |. we here used a grid 
with 128^ points. In one set of simulations, the situa- 
tion was idealized, by using periodic boundary conditions 
and perfectly identical master and slave systems. In the 
other set, circular boundaries were used. Moreover, we 
tried to include typical experimental imperfections, such 
as speckles and spatial variations of the nonlinear sen- 
sitivity. Since it is very difficult to quantify the experi- 
mental imperfections exactly, these simulations give only 
a qualitative picture. To simulate speckles, we added 
white spatial noise to the pump beam (2% of the total 
intensity in non-zero Fourier modes). The nonlinear sen- 
sitivity was given smooth spatial profiles, different for 
master and slave, with an variation of ±1% (RMS) 20]. 

The resulting cross-correlation coefficients are shown in 
Fig. |S1 Even though the pump was set to 7 times thresh- 
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FIG. 5: Cross-correlation coefficient cab from numerical sim- 
ulations, for idealized systems at 7 times threshold (O), and 
with imperfections included at 4 times threshold (a). 

old for the idealized case, the correlation grows quite 
fast with the coupling strength and reaches almost unity. 
Whereas when imperfections are included, the synchro- 
nization is still present, but the cross-correlation stays 



significantly smaller. 

In conclusion, we have found a first experimental evi- 
dence for the synchronization of spatiotemporal chaos in 
a nonlinear-optical system, realized in a unidirectionally 
coupled master-slave configuration. The dependence of 
the cross-correlation between master and slave state on 
the coupling strength clearly indicates synchronization. 
Corresponding numerical simulations reveal how spatial 
inhomogcneities counteract synchronization. 

The underlying mechanisms of this synchorization 
(from x'-^pi-oi'- place) of spatiotemporal disorder appear 
to be quite general; an observation in other nonlinear ex- 
tended systems should be possible. Good candidates, be- 
sides other nonlinear optical systems, are photosensitive 
chemical reaction-diffusion systems, wich can both be de- 
tected and manipulated with light waves. A necessary 
translation between detection and controlling light or an 
amplification of the detection light can be performed by 
means of optically addressable spatial light modulators 
[2ll |. The present work can also be considered as a step 
to secure parallel communication, as proposed in . 
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